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We study a single self avoiding hydrophilic hydrophobic polymer chain, 
through Monte Carlo lattice simulations. The affinity of monomer i for water 
is characterized by a (scalar) charge Aj, and the monomer- water interaction 
is short-ranged. Assuming incompressibility yields an effective short ranged 
interaction between monomer pairs (i, j), proportional to (Aj + Xj). In this 
article, we take \ = +1 (resp. (A, = — I)) for hydrophilic (resp. hydropho- 
bic) monomers and consider a chain with (i) an equal number of hydro-philic 
and -phobic monomers (ii) a periodic distribution of the Aj along the chain, 
with periodicity 2p. The simulations are done for various chain lengths N, 
in d = 2 (square lattice) and d = 3 (cubic lattice). There is a critical value 
p c (d, N) of the periodicity, which distinguishes between different low temper- 
ature structures. For p > p c , the ground state corresponds to a macroscopic 
phase separation between a dense hydrophobic core and hydrophilic loops. 
For p < p c (but not too small), one gets a microscopic (finite scale) phase 
separation, and the ground state corresponds to a chain or network of hydro- 
phobic droplets, coated by hydrophilic monomers. We restrict our study to 
two extreme cases, p ~ O(N) and p ~ 0(1) to illustrate the physics of the 
various phase transitions. A tentative variational approach is also presented. 
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I. INTRODUCTION 



A very popular approach to the protein folding problem, is to emphasize the heterogeneity 
of a protein due to the different side chains |T[ . Of special importance in this context is the 
model of a (quenched) random hydrophilic hydrophobic chain, since (i) proteins are usually 
designed to work in water (ii) the first step of the folding transition may correspond to the 
formation of a hydrophobic core H. In this model, each monomer is described by a single 
"charge" Aj, and the polymer- water interactions are modelled through the Hamiltonian: 

N Af 

u pw = - X) _ Ra ) ( 1 ) 

i=l a=l 

where a(f) denotes a short-ranged (Van der Waals-like) monomer- water molecule interaction, 
iV and M denote respectively the number of monomers and of water molecules, and fj and 
R a , their respective positions. Assuming that the system is incompressible, one gets JUII 

i N N N 

n pw = +-ED^ + \) a ^ - ?j) ~ A Y, x i ( 2 ) 
^ i=i j=i i=i 

where A = J2r a {^)- The second term in (fj) is a constant for the quenched disordered 
chain (as well as for the periodic chain studied in this article). It will therefore be omitted 
henceforth. The phase diagram of this random hydrophilic hydrophobic polymer has been 
studied by various methods (mean-field, dynamics, replica variational calculations,....) for 
simple disorder distributions of the Aj 's The results can be roughly summarized 

as follows 

(i) for strongly (on average) hydrophobic chains, one expects a 9 collapse transition to 
occur first, followed at lower temperatures, by a scale dependent freezing transition. 

(ii) for chains that are (on average), either weakly hydrophobic or hydrophilic, one ex- 
pects a freezing transition that is perhaps characterized by a "random first order" transition 

§• 

An interesting application of this model to proteins has recently been published 0, even 
though their hydrophobic content appears not to be random [ fTU| . 
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As a first step towards the protein folding problem, we will study numerically a single 
self avoiding chain, where the distribution of the Aj 's is periodic along the chain. To be 
more specific, the chain of iV monomers is made out of periodically alternating blocks of 
p hydrophilic monomers (A, = +1) and p hydrophobic monomers (Aj = —1). Following 
equation (Q), the Hamiltonian of the chain is defined as 

i N N 

« = +2 + Aft-f.,-) (3) 

i=l j=l 

where A(f*j — fj) denotes a lattice A function: A(r) = 1, if r connects nearest neighbours, 
and otherwise. In this paper, we will consider square (d = 2) and cubic (d = 3) lattices. 

The parameter p is an important ingredient of the problem since there is a critical 
value of p, that we note p c (d,N), above which the equilibrium ground state consists of 
a single hydrophobic core. The argument is as follows; for a spherical core (d = 3), the 

2 

number of points on the surface is of order iVa, whereas the surrounding hydrophilic loops 
occupy a number of surface points of order ^. Therefore, one expects a single macroscopic 
hydrophobic core for p > O(N^), and many microscopic hydrophobic droplets for p < 
O(N^). We therefore get p c (3,N) ~ O(N^). For d = 2, the steric constraints are much 
stronger, and one hydrophilic loop anchored on the hydrophobic core screens part of its 
"surface". Since the distance between the two ends of the hydrophilic loop typically scales 
with the radius of gyration of the core, that is also with its "surface" , this screening effect 
is very strong. It implies that p c (2,N) ~ O(N). In the case p < p c (d,N), corresponding 
to many hydrophobic droplets, one should further distinguish between a linear or branched 
topology of the hydrophobic droplets. For very small p, the very formation of these droplets 
is impeded by the hydrophilic parts of the chain (see below) 

The numerical simulations for the periodic chain of equation (|J) are done on the square 
or cubic lattice, using the Multiple Markov chain method The plan of the paper is 

the following. We first recall, in section [TT] the principles of the Multiple Markov chain 
method in Monte Carlo simulations of polymer chains. As a preliminary test, we apply this 
method to the two dimensional 9 collapse of a purely hydrophobic chain. We then consider 
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the periodic hydrophilic hydrophobic chain in two extreme cases, namely p ~ O(N) and 
p ~ 0(1), to capture the physics of the various phase transitions. Putting numbers on the 
above estimates of p c shows that it will be numerically convenient to study the single core 
phase in d = 3 and the multiple cores' phase in d = 2. In section |TTT| , the case p ~ 0(N) is 
studied in detail for d = 3; we also briefly consider the corresponding transition in d — 2. In 



section |TV|, we study the case p ~ 0(1) which, as stated above, has a multiple droplet low 
temperature structure: for d = 2, we present evidence for the existence of an intermediate 
branched phase, if p is not too small. The d = 3 case requires very long chains and is only 
briefly considered. Finally, we also present a tentative variational method for this case. 

II. THE MULTIPLE MARKOV CHAIN METHOD 

A. Summary of the method 

In this section we give a quick description of the numerical techniques which we use to 
calculate thermodynamic (and/or geometric) properties of the chains as a function of the 



temperature. These techniques have been discussed in detail in reference [12], and can be 
summarized as follows. 

The implementation of the Metropolis Monte Carlo method relies on the multiple Markov 
chain sampling. First, one generates with the simple Metropolis heat bath, a Markov chain 



at temperature T: this procedure makes use of an hybrid algorithm, which has pivot [13 



as well as crankshaft moves fl4]]. Pivot moves are of a global type, and operate well in 



the swollen phase; crankshaft moves, which are of a local type, are essential in speeding up 



convergence close to the collapse phase transitions fl2]l . In these calculations, each Monte 
Carlo step consists of 0(1) pivot move and 0(N) crankshaft moves. 

One may then run in parallel a number m (typically 15 — 20) of these Markov chains at 
different temperatures T x > T 2 > . . . > T m . We allow the chains to interact (by possibly 
swapping conformations) as follows. Among the m chains, we select, with uniform proba- 
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bility, two chains (a,/?) at respective temperature T a and Tp. A trial move is an attempt to 
swap the two current conformations of these chains. If hk{T) is the probability of getting 
state K at temperature T (that is ttk(T) ~ e ~, in obvious notations), and S a and Sp 
are the current states in the a-th and (3-th chain, then we accept the trial move (i.e. we 
swap S a and Sp) with probability 

r(S a ,Sp) = mm [l p \. (4) 

Note that the whole process is itself a (composite) Markov chain. Since the underlying 
Markov chains are ergodic, so is the composite Markov chain. Furthermore, the composite 
chain obeys detailed balance, since the swap moves as well as the moves in the underlying 
chains obey detailed balance [0. The swapping procedure dramatically decreases the cor- 
relations within each Markov chain, and produces little CPU waste since, in any case, one 



is interested in obtaining data at many temperatures [I I 



B. The two dimensional 6 transition 



We first consider the application of the Multiple Markov chain method to the two dimen- 



sional 6 collapse transition of a purely hydrophobic chain (see ref [^2| for the case d = 3). 
We are mostly interested here in locating the thermodynamic 8 transition. 



We show in Figure 1 the plot of the average scaled radius of gyration (4tL^) versus 



temperature, for different values of N. Finite size scaling theory [|ll^, |il||T7j ,|T 
that 



<R 2 > 



f{tN*) 



predicts 



(5) 



where t 



I T — Tg 

1 T e 



is the reduced temperature for N large (T e = T e (N = oo)), $ a crossover 



exponent, and f(x) a scaling function with a finite value for x — > 0. Using the exact result 
2uq = | yields an estimate of the 9 temperature: Tq ~ 1.5, in agreement with recent 
estimates To study the critical behaviour of the specific heat close to the 6 point is a 



6 



notoriously difficult problem in d = 2. We have also verified this point, and we have only 
extracted from the N dependence of the peak of the specific heat (Figure 2) a large N critical 
temperature T' e ~ 1.5 in broad agreement with the value obtained from the behaviour of the 
radius of gyration. 



III. THE CASE P =0(N) 

A. Numerical simulations for d=3 

This case corresponds to a single hydrophobic core ground state since p c (3, N) ~ O(Ns). 
In the case of the 9 collapse transition, the iV monomers of the chain play an equivalent role 
(neglecting end effects for a long enough chain). In the present problem, the repulsive inter- 
action between hydrophilic monomers leads us to consider two possibilities (i) the collapse 
transition first occurs in a single hydrophobic block of length p (ii) the collapse is due to a 
cooperative effect of the ^ hydrophobic blocks, in a way similar to the 6 transition. Scenario 
(ii) is consistent with a unique phase transition of a discontinuous character occurring at 
temperature T^y). Scenario (i) is a priori consistent with a collapse transition at Tg(p), and 
raises the question of a (surface induced) sticking transition of the individual hydrophobic 
blocks. 

In either case, the collapse transition is expected to be discontinuous, with a jump of the 
radius of gyration. We have fixed p = y in the simulations, and let iV vary from 80 to 640. 
A typical low temperature configuration for iV = 640 is shown in Figure 3: as expected, it 
displays a macroscopic phase separation between the hydrophobic and hydrophilic parts of 
the chain. To get an estimate of the critical temperature, we have first considered the radius 
of gyration < R 2 > of the complete chain. As shown in Figure 4, the exponent of < R 2 > 



is, at all temperatures, given by the self avoiding walk (SAW) value (vsaw — 0.588 [p^ ,p3f| ) . 

Our data for the large values of N are in agreement with scenario (i), that is a single 
hydrophobic block collapse, since the critical temperature T c is very close to Tg(p) (see Table 
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I). A better estimate of T c comes from specific heat measurements, since the specific heat 
has a quite spiky character (Figure 5), unlike its 9 point equivalent [jT^] . We again get an 
estimate very close to T${p). We have also tried a finite size scaling analysis to find T c , 



for large N (and therefore large p). Following a well established path (see e.g. |M| and 
references therein), we have considered the N dependence of the height and position of the 
peak in the specific heat. At the transition, and in the scaling limit, one expects the height 
peak C* to scale like 

C* ~ iV 2 *" 1 (6) 



where the crossover exponent $ and the specific heat exponent a are related |25|] through 
the relation $(2 — a) = 1. At the three dimensional 9 transition, one has $ = 0.5 and a = 0. 
At a (thermal) first order transition, one has a — 1, yielding C* ~ N. Finite size scaling 
also implies a critical temperature shift AT C = T c (oo) — T C (N) ~ -^L-, yielding AT C ~ 
for a first order transition. Figures 6(a) and 6(b) show our results, and confirm that our 
simulations are not done in the scaling regime. It is well known indeed that, in order to get 
a precise estimate of the thermodynamic 9 temperature, one has to study very long chains 
(typically iV > 1500). For our problem, since p = y, we should study the case N = 12000, 
which is presently out of reach. We have done an independent simulation with N = 800, and 
p = ^ = 200. The peak in the specific heat occurs for T c ~ 2.5, very close with Tg(p = 200) 
(see Table I). 

We therefore believe that, in the thermodynamic limit (N,p — > oo), one has a discon- 
tinuous single block collapse transition; the critical temperature T c (oo) is the same as the 
collapse temperature of the fully hydrophobic chain (Tg = T 9 (oo) ~ 3.7 [pH|). The transition 
is well characterized by the phase separation order parameter |6].[7| : 

< 5R 2 >=< R 2 phil > - < R 2 phob > (7) 

where < I& ha > (resp. < R 2 phob > is the squared radius of gyration of the hydrophilic 
(resp. hydrophobic) monomers. In Figure 7, we have plotted the order parameter ( ^iJf^ ) 
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as a function of temperature, for various values of N. Its behaviour is consistent with the 
previous results. 

An intriguing feature of the single block collapse mechanism is that it seems to imply the 
existence of an intermediate phase, between the low temperature phase depicted in Figure 
3, and the swollen coil phase. This intermediate phase is a necklace (or network) of single 
hydrophobic blocks, and its free energy differs from the low temperature free energy by 
a surface free energy. We have not found this intermediate phase in our multiple chain 
Monte Carlo Method, and this may have several causes. It is for instance possible that the 
temperatures of the various chains T\ > T2 > . . . > T m of our simulations had too large a 
spacing to find this phase. Further work, in particular a precise determination of surface 
properties, is needed on this point. 



B. Numerical simulations for d=2 



In this case, one expects p c (2,N) ~ O(N), yielding two different situations. If p > p c , 
one will get a single hydrophobic core below the collapse transition (Figure 8). One may 
then study the (full) radius of gyration at high and low temperatures (Figure 9), obtaining 
in both cases self avoiding behaviour (i.e. v = | |26fl). If p < p c , several hydrophobic 
cores appear (Figure 10). On the square lattice, we have found numerically that 0.06 < 
p c (2, N)/N < 0.08. One may get a feeling for this result by considering a square shaped 
hydrophobic single core that is fully surrounded by hydrophilic loops: the above estimates 
corresponds to a total number of hydrophilic loops approximately equal to 6 — 8. This 
argument is only suggestive since the single core becomes elongated as p — > p+. This is 
clearly due to the screening effect of the hydrophilic loops: the hydrophobic core tries to 
maximize its perimeter at fixed surface (Figure 11). As for the phase transition for p > p c , 
one may say that specific heat data display a rather smooth behaviour; the transition is not 
markedly discontinuous. As already mentionned, extracting a more detailed information 
from these data is rather tricky in two dimensions. 
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IV. THE MANY HYDROPHOBIC DROPLET CHAIN P=0(1) 



A. Numerical simulations for d=2 

The large value of p c (2, N) shows that the multiple droplets' phase should be a priori 
easier to study in two dimensions. We show in Figure 12 a typical low temperature configu- 
ration, which displays branched polymer features |[27||. From section [111 B| we know that the 



maximum number of hydrophobic monomers in a droplet is Umax ~ 12 — 16 p. To further 
investigate the branched character of the phase, we have studied the case p = 8, with N 
ranging from 80 to 1200 (other values of p are briefly considered below). For p = 8, the 
high and low temperature exponents v of the radius of gyration are shown on Figure 13. 
Above the transition, we get SAW behaviour; below the transition, we get v ~ 0.64, which 



is indeed close to the branched polymer (BP) value [28 . 

It turns out that low temperatures are difficult to study because of non-equilibrium 
effects, so that we are not able to follow in detail the thermal evolution of the branched 
phase. This is partially due to the fact that the Monte Carlo method of section |TJ has not 
been optimized to deal with branched phases. Another reason may be the possible existence 
of a dynamical phase transition towards some kind of glassy branched state (see section 
IV □ ). We have therefore restricted our study to the phase transition between the high 



temperature (SAW) and low temperature (BP) phases. Using the same finite size argument 
as in section |11B| , we plot in Figure 14 the scaled radius of gyration (4^^) vs temperature. 
Various values of the unknown exponent v* have been considered (Figure 15). Our results 
show the existence of a phase transition at T c ~ 0.8 ± 0.1, and evidence for a new critical 
behaviour (v* ~ 0.70 ± 0.03) at T c . The phase transition seems also to be quite smooth, 
if one considers specific heat data. Another "experimental" observation concerns the size 
distribution of the hydrophobic droplets: below the transition, we find that most of the 
droplets do not reach the maximum size allowed n M Ax- This can be interpreted as an 
entropic effect, very much along the lines of reference ||. 



10 



We conclude this section by a few remarks on the role of p. We have also considered 
the case p — 4, and p = 10, 12, with the same range of values of N. We do not find a 
clear evidence for a phase transition for p — 4, whereas we find evidence for two transitions 
for p — 10, 12. In the latter case, the branched phase gives way at low temperature to a 
reentrant self avoiding chain of finite hydrophobic droplets. This shows that the balance 
between linear and branched topologies is very dependent on the value of p. If p is too 
small, the formation of the droplets is impeded by the repulsive hydrophilic monomers. For 
large p, a local collapse is possible, favoring the linear topology at low temperature. For 
intermediate values of p, a non local hydrophobic collapse is apparently favored, yielding a 
branched topology. These issues will be further tackled in section [IV C| . 



B. Numerical simulations for d=3 

A typical low temperature configuration is shown in Figure 16, for p = 4, N = 720. The 
properties of the multiple hydrophobic cores phase are difficult to study, since one needs very 
large values of N. Furthermore, for d = 3, a SAW at the 9 point and a branched polymer 
(BP) have the same exponent u$ = v^p = 0.5 [25|, which makes a detailed scaling analysis 
difficult. 



C. Variational method for the many droplet phase 

Following traditional polymer physics methods [|IJ, we will study the low temperature 
branched phase in a variational way. The basic steps can be summarized as follows: 

(i) one derives an effective quantum Hamiltonian. 

(ii) ones uses a ground state approximation, together with a saddle point approximation. 

(iii) finally, one performs a variational calculation, and minimizes the free energy with 
respect to the relevant parameters. 

Steps (i) and (ii) are familiar in the context of the usual 9 collapse transition. Since one 
is then interested in a macroscopic collapse, a continuous description of the chain is valid; 
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furthermore, the ground state approximation holds for long enough chains, since there is a 
bound state representative of the collapsed globule. Finally, step (iii) is usually implemented 
with a constant or Gaussian density around the center of mass of the collapsed globule, which 
can be taken as fixed in all calculations. 

On the contrary, what we have in the low temperature branched phase is a inhomogeneous 
liquid of microscopic hydrophobic droplets. If one follows the above procedure, one must 
take the extensive entropy of these droplets (i.e. the degeneracy of the saddle point) into 
account. Since we believe that our approach may be of interest in other contexts f|, we will 
assume that a continuum description of the chain is valid, and derive the simplest form of 
the associated Hamiltonian. We will also assume that ground state dominance holds. 

The partition function of the hydrophilic hydrophobic chain reads 



d 



Y[dn e-i^E^+i-^) 2 -^ ( 8 ) 

i 

where a is a typical monomer length and i = 1,2,...N. The lattice Hamiltonian (7i) has 
been derived in equation ([$]), and its off-lattice version reads 

FH = \ E K + P(\ + A,)] Sin - r 3 ) 

+ 1 £ S{fi - fj) 5(fj - f fc ) 

1 

+^7 £ y° 5 (fi - fj) 5(fj - r fc ) s(f k - fi) (9) 

Note that we have also included three (u> ) and four (y ) body terms for reasons that will 
soon become clear. We also assume (fi = ?n = 0). 
Defining the local density p(r) as 

p(f) = Y, sf-fi) (io) 

i 

we have 



Z = J V(f){f)Vp(r) C{p, (p) exp (i J d d r<p(r)p(r) - J d d r 



(ii) 
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with 

i 

For the periodic hydrophilic hydrophobic chain, we introduce transfer matrices T± for Aj = 
±1 and get 

C(p,0)=<O|[T^]i|O> (13) 



with 



and 



< f|T_|f >= e -5^) a -^0+/W> (15) 

In equation (|13|), the differences between p ~ O(N) andp ~ 0(1) are clearly displayed. From 
now on, we will set p = 1, without questioning any further the existence of the continuum 
limit in this case. Using the identity 

< r\T + T„\r >= J d d r < r\T + \r >< f \T + \r > (16) 

together with equations (|14|) (|15|) leads, to lowest non trivial order in <j)(f) and p(r), to 

< f\T + T^ >= e -^- r ^ 2 - 2l ^ + ^ 2 ^p) 2 (r) (17) 



which implies 



N 

C(p,0) = Tr exp-(-Ho) (18) 



with 



a 2 - a 2 



n ° = ~4d v + 2%m ~~ Td p{yp){f) (19) 

It is quite clear that our derivation is not rigorous. We nevertheless feel that the p{f) de- 
pendent term of the Hamiltonian 7i is physically sound since it favors inhomogeneous high 
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density regions (droplets). A better approximation would presumably involve higher deriva- 
tives of p(r), which clearly define typical droplet sizes. Assuming ground state dominance 
in (|T9|), and performing a saddle point approximation on <ft(r) in (JTTJ) yields 

p(f) = N^ 2 (f) (20) 

where \l/(r) is a normalized wave function. We then obtain the "free energy" per monomer 
as 

^ = mm / d d r G(*(r)) (21) 

where 

+i J d"r*(rl (-^V 2 - ^(V(Nr-)m^j 9(f) (22) 

At this point, it is useful to remark that one needs in this approximation to introduce four 
body interactions, as in the disordered case, and for the same reasons H. The fact that 
attractive multi-body interactions in homopolymers may induce a SAW-BP phase transition 
has been previously noted for a specific model in reference [30]. 



At low temperature, the last term of equation ([220 plays a dominant role. This term, 
as mentionned above, tends to create an interface between hydrophilic and hydrophobic 
regions. If one uses a variational wave function *f?o(f) given (for d — 2) by 

%(r) = a cos(— ) cos(-^) (23) 

the normalization condition implies 

a 2 R 2 ~ 0(1) (24) 

where R is the linear dimension of the system. Plugging this estimate in equation fl2~2|) shows 
that a low free energy is obtained for I ~ a and a finite average density p av = ~ 0(1). In 
other words, the low temperature phase obtained from our variational approach is a dense 
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phase (since v = |, with d = 2), made of microscopically phase separated regions. In 
section ( [IV A| ) , we obtained from our simulations the value v ~ 0.64 below the the (SAW)- 
(BP) phase transition. This result is not compatible with the value v = 0.5 that we get 
through the variational method. Some possible explanations for this "discrepancy" are as 
follows 

(i) The simulations were done for p = 8 on a lattice, and the variational method was 
applied to the case p = 1, within a continuum limit approach. 

(ii) A major difference between the present multi droplet collapse and the 9 single core 
collapse is that one has to take into account the degeneracy of the saddle point equation for 
<f)(r) in evaluating the "true" free energy. In other words, there is a droplet entropy that 
must be considered. Sticking with the variational function of equation (|23|), it is easily seen 
that this entropy favors large values of -. A precise calculation is difficult, and we will not 
comment upon this point anymore. 

(iii) We believe however that the main reason for the difference between the exponents 
v stems from the use of a ground state approximation in estimating the right handside of 
equation flT8p. As far as we know, this approximation, which relies on the existence of a 
bound state in the Hamiltonian Tio, works well for dense (finite density) phases. It does not 
a priori describe a branched polymer, which has a vanishing density. Physically, a dense 
phase is not favorable because of the repulsion between the hydrophilic monomers. 

Altogether, our results seems to indicate that the finite p chain may undergo zero (p = 4), 
one (p = 8) or two (p = 10, 12) phase transitions. Ground state dominance is never a 
valid approximation, since one deals with either linear [y = 0.75) or branched [y ~ 0.64) 
structures. Reentrant behaviour, similar to the one described for p = 10, 12, has been 



found recently in related models pl| , p2| . Another issue of interest concerns a possible low 
temperature dynamical (glass) transition, similar to the one described in reference M. 
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V. CONCLUSION 



We have studied a periodic hydrophilic hydrophobic chain. An important ingredient 
of the physics of this problem is the value of the period 2p. The low temperature phase 
consists of a single (p > p c (d,N)) or multiple (p < p c (d,N)) hydrophobic core(s), where 
p c (3, N) ~ O(N^) and p c (2,N) ~ O(N). Using Monte Carlo calculations, we have studied 
the case p ~ O(N) and indeed found a macroscopic phase separation between the two types 
of monomers in d — 3, and two possible regimes in d = 2. The second case, (p ~ 0(1)) 
yields for both dimensions a low temperature phase, consisting of a chain or network of 
microscopic hydrophobic droplets linked by hydrophilic filaments. We have studied this 
phase numerically in d = 2. A connection with both the branched polymer chain [27] and 
the random hydrophilic hydrophobic chain H,|5|j6|j7|,pl[1 is physically appealing and shows up 
in the tentative variational treatment of this phase, as given in section [IV C| . As for proteins, 



the existence of periodic hydrophobicity patterns in secondary structures [TO]] suggests that 
our model may have some relevance in explaining the typical size of single domain proteins 
(N ~ 120 — 150 residues). Further work in these directions is in progress. 

It is a pleasure to thank Henri Orland for fruitful discussions and suggestions, and 
Bernard Derrida for interesting comments. 
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Figure captions 

Figure 1: Scaled radius of gyration ( ) vs temperature for the purely hydrophobic 
chain in d = 2, for N = 80 (A), 160 (D),240 (O),480 (x),640 (+). A crossing occurs for 
T e ~ 1.5. 

Figure 2: Specific heat vs temperature, for the purely hydrophobic chain in d = 2, for 
the same values of N. The extrapolated critical temperature is T' B ~ 1.5. 

Figure 3: Typical phase separated configuration (d = 3, p = ^, N = 640). Circles 
denote hydrophobic monomers. 

Figure 4: Log-Log plot of < R 2 > vs iV at various temperatures for (d — 3, p — y), 
and JV = 80, 160, 240, 360, 480, 640, 800. The temperature is T = oo (□), 3.33 (A), 2.0 (*), 

,1.0 (O),0.5 (+). The upper and lower straight lines have slopes compatible with self 
avoiding behaviour 2p S aw — 1-176. 

Figure 5: Specific heat vs temperature for (d — 3, p — y), and iV = 80 (O), 160 (□), 

240 (O),360 (x),480 (+),640 (A). Note the increase of the peak as well as its shape, 
when N increases. 

Figure 6: (a) Plot of the specific heat peak C* vs iV for (d = 3, p = ^), and iV = 
80, 160, 240, 360, 480, 640 (b) Plot of the critical temperature vs for the same parameters. 
Error bars correspond to one standard deviation. 

Figure 7: Scaled phase separation parameter (see equation (^)) vs temperature, for 
(d = 3, p = f ), and iV = 80, 160, 240, 360, 480, 640. 

Figure 8: Low temperature phase separated configuration (d — 2, p — ^, N — 240). 
Black triangles denote hydrophilic monomers. Note the isotropic shape of the hydrophobic 
core (p » p c ). 

Figure 9: Log-Log plot of < R 2 > vs iV at various temperatures for (d — 2, p — y), and 
N = 80, 160, 240, 360, 480, 640, 800. The temperature is T = oo (□), 1.4 (A), 1.0 (*), 0.5 (O). 
The upper and lower straight lines have slopes compatible with self avoiding behaviour 
Zvsaw = 1-5. 



19 



Figure 10: Low temperature multiple cores' configuration displaying the screening effect 
of the hydrophilic loops (d = 2, p = g, TV = 240). 

Figure 11: Low temperature phase separated configuration (d — 2, p — || ~ p c +, 
N = 240). Note the elongated shape of the hydrophobic core. 

Figure 12: Typical multiple cores' configuration (d — 2, p — 8, N — 1200). Black 
triangles denote hydrophilic monomers. 

Figure 13: Log-Log plot of < R 2 > vs N at various temperatures for (d — 2, p — 8), and 
N = 80, 160, 240, 360, 480, 640, 800, 1200. The temperature is T = oo (□), 1.0 (A), 0.7 (*), 

,0.5 (O),0.33 (+).The upper straight line has a slope 2usaw = 1-5, whereas the middle 
line has slope 2u B p, where v B p — 0.64 is the branched polymer (BP) value. The bottom 
line corresponds to the collapsed value 2u — 1. 

Figure 14: Scaled radius of gyration ( ^1? ) vs temperature for (d — 2, p — 8) and 
iV = 80 (D),160 (O),240 (x),480 (+),640 (A). A crossing occurs for v* ~ 0.70, yielding 
T c ~ 0.8. 

Figure 15: With the same data, no crossing occurs for (a) v* = i>saw — 0.75 (b) 
v* = v e = ± (c) v* = v BP ~ 0.64. 

Figure 16: Typical multiple cores' configuration (d = 3, p = A, N = 720). Circles 
denote hydrophobic monomers. 
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Table caption 



Table I: Comparison between the critical temperature of the hydrophilic hydrophobic 
chain T c = T c (p, N), and the 6 transition temperatures [|T^] of a fully hydrophobic chain of (i) 
p monomers (ii) y monomers. The first three lines correspond to p = y (N = 80, 640, 800). 
The last line corresponds to p = ^ (N = 800). The value Tg(p = 10) has been obtained 
using exact enumeration techniques. 
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